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ABSTRACT 

We present the results of a Doppler tomographic reconstruction of the UV 
spectra of the double-lined, O-binary HD 100213 based on observations made 
with the International Ultraviolet Explorer (IUE). We used cross-correlation 
methods to obtain radial velocities, confirm the orbital elements, estimate the 
UV flux ratio, and determine projected rotational velocities. The individual com- 
ponent spectra are classified as 07 V + 08 V using UV criteria defined by Penny, 
Gies, & Bagnuolo. We present a model fit of the eclipsing light curve from obser- 
vations from the HIPPARCOS satellite and published observations of Andersen 
& Gronbech. We derive an orbital inclination, i = 77°.7 ± 1?0. This analysis in- 
dicates that both stars are currently experiencing Roche lobe overflow (RLOF), 
which confirms earlier results that this is one of only a few massive contact bi- 
naries. Our derived masses, M p /M Q = 16.8 ± 0.4 and M S /M Q = 10.5 ± 0.3, are 
significantly lower than those computed from the Doppler shifts of lines in the 
optical spectrum. We suggest that the difference occurs because mutual irradia- 
tion decreases the upper atmospheric temperature gradient in the inwards facing 
hemispheres of both stars, which makes lower excitation lines appear weaker ther 
and shifts their center-of-light away from the center-of-mass. We compare the 
current state of HD 100213 with predicted outcomes of massive close binary evo- 
lutionary models, and suggest that the system is currently in a very slow Case 
AA mass transfer stage. 

Subject headings: Stars: binaries: spectroscopic — stars: binaries: eclipsing — stars: 
early-type — stars: fundamental parameters — stars: individual (HD 100213) - 
ultraviolet: stars 
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Introduction 



In this series of papers we have utilized the tomography algorithm to better study 



(Baenuolo & Gies 


1991 


), d> Per ( 


Thai 


ler et al. 


1995 


), HD 15224* 


1 ( 


Pennv. Gies. & Baenuolo 


1999 


), HD 135240 


(Pennv et al. 


2001 


), and 5 Ori ( 


Harvin et al. 


2002 


). This time we turn 



our atte ntion to the do uble-lined, eclipsing, massive binary syst em, HD 100213 (TU M us; 



08.5 V, IWalborn 



(11975 



19711 ) . This binary was originally studied by I Andersen fc Gronbech 



AG 75) who presented a combined spect roscopic and ph o tome tric solution and 



Stickland et al. 



(I1995I . ST95) performed a 



concluded that the stars were in contact. Later 
spectroscopic study using 24 high-resolution UV spectra from the International Ultraviolet 
Explorer (IUE). Although their orbital period agreed well with that from AG75, their orbital 
semiam plitudes (K) were s ignificantly smaller, resulting in smaller binary separation and 



masses. 



Terrell et al. 



( 2003 



T03), with an eye towards resolving this discrepancy, obtained 
new optical (AA3900 — 9200 A) spectra at two orbital phases close to quadrature, as well 
as new photometric observations. These they combined with the earlier spectroscopic and 
photometric data of AG75 and HIPPARCOS observations. The results of this combined 
analysis supported the AG75 solution and appeared to rule out the smaller amplitude 
values from ST95. As both ST95 and T03 note, that this discrepancy between UV and 
optically determined K values is not unique to TU MUS, but also is seen in several other 
massive semi-detached and contact binaries. 



Linder et al 



(120071 . L07) obtained optical spectra of TU Mus to determine whether the 
system displayed the Struve-Sahade (SS) effect, i.e., a change in the relative line depths 
of the components between quadrature phases. While they did not observe this effect, 
they did note several issues that are illuminating. First, their orbital elements (specifically 
semiamplitudes) vary depending upon which lines were used, i.e., the K values from He II 
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lines were smaller than those from He I. Second, the equivalent widths of two He I lines 
are phase locked with the orbital motion. They are at maximum strength when the back 
(opposite the companion) of the star is facing towards us and at a minimum half a orbit 
later. They conclude that the He I lines preferentially form over a portion of the stellar 
surface, and possibly an analogous, but opposite, effect may be occuring with the UV lines, 
due to mutual irradiation of the inner sides of the stars. 

As one of only a few contact massive binaries, the TU Mus system is an important 
test case for interactive binary evolution models. It is critical that well-determined values 
for the current separation and component masses be determined. In this paper we present 
our own analysis of the IUE observations of HD 100213 resulting in a double-lined orbital 
solution (§3), our tomographic reconstruction of the composite spectra into separate 
primary and secondary spectra and their respective spectral classifications (§5), individual 
projected rotational velocities (§4), and our model of the eclipsing light curve (§6). We 
discuss the discrepancy between the orbital semiamplitudes determined from optical spectra 
versus those from UV spectra in §7. We compare the individual masses based on the 
spectroscopic and photometric orbits to those predicted by theoretical methods and discuss 
their implications (§8). 



2. Observations and Reductions 

There are 24 high dispersion, short wavelength prime camera spectra of HD 100213 
available from IUE. These spectra were obtained from Multi-Mission Archive at Space 
Telescope Science Institute {MAST). One spectrum, SWP 54386 was found to be flawed, 
with several large gaps in spectral coverage. ST95 note that this spectrum was severely 
overexposed due to a satellite tracking problem. This spectrum was not used for our 
analysis. The individual SWP image numbers and heliocentric Julian dates of mid-exposure 
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of the remaining 23 observations are presented in Table 1. 



EDITOR: PLACE TABLE [T] HERE. 



The spectra were manipulated in several stages to produce a matrix of spectra (in 
dimensions of wavelength and time) rectified using a common set of relatively line-free 
zones and sampled with a uniform log A wavelength grid. The major interstellar absorption 



lines were replaced by st r aight line segments in the processing. Details are given in 
Penny. Gies. fc Bagnuolol ( 1997 ). 



3. Radial Velocities and Orbital Elements 



Radial velocities and orbital elements based on this same set of spectra were presented 
by ST95. However, we decided to check their results using our suite of cross-correlation 
techniques to obtain rad ial velocities and e stimates of projected rotational velocities and 



flux ratio (for details see 



Penny et al 



19971 ). We used two spectra (SWP 22108, SWP3 7429) 



Penny 



19961 ) for 



of HD 34078 (09.5V; projected rotational velocity Vsini = 26.5 km s _1 
cross-correlation with the spectra of HD 100213. The relative velocities were determined 
by fitting two Gaussians to the composite cross-correlation functions (ccfs). The widths of 
the two Gaussians and their intensity ratio were determined by least-squares fits of ccfs 
with well-separated peaks, i.e., those obtained within ±0.12 of the quadrature phases. The 



final radial velocities (Table 1) were transfor med 



radial velocity of HD 34078 (+54.4 km s" 1 , 



Gies 



rom r el ative to absolute by adding the 



19871 ). 



Qian et al. 



(120071 ) investigated 



possible cyclical variations in orbital period of TU Mus and suggested that there is a 
tertiary companion. The predicted low mass of this star, « 2.OM , indicates that it will 
not significantly contribute to the IV E spectra. In addition there is a well known faint 
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companion, with a magnitude V = 13.4, that is located 15" to the southwest of TU Mus, 
but we see no evidence of any tertiary spectrum in our ccfs. 



We determined the orbital elements using the program (SBCM) of 



Morbey fc Brosterhus 



(119741 ) . This program is limited to individual solutions for the component stars. We 
assigned zero weight to those velocities obtained from spectra within about 0.15 phase of 
the minima of eclipse and 1.0 weight on all other observations. The orbital elements based 
on the IUE measurements are presented in Table 2. The combined primary and secondary 
solutions of ST95 and T03 are presented along with the separate primary and secondary 
solutions from this paper. We note that T03 performed a simultaneous solution with both 
the spectroscopic and photometric data using the 2003 Wilson-Devinney (WD) code. They 
do not list their resulting K values, but do present their final mass ratio, period, inclination, 
and separation, from which we determined their individual semiamplitudes (shown in Table 
2). The numbers in parentheses refer to the errors in the last digit quoted. We ran the 
program holding the eccentricity at zero, and allowed the other parameters to vary. We 
reached excellent agreement with ST95 in all parameters. The radial velocity curves based 
on the data are illustrated in Figure 1. The velocity residuals and orbital phases from the 
individual fits are listed in Table 1. The epoch To corresponding to zero phase corresponds 
to the radial velocity maximum of the primary star. 

EDITOR: PLACE TABLE HERE. 



EDITOR: PLACE FIGURE Q] HERE. 



Our results are not in agreement with those from T03 who find significantly higher 
orbital semi-amplitude velocities for both stars. As noted above the disagreement between 
orbital semi-amplitude values determined from UV spectra and those from optical spectra 
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is not unique to TU Mus. Four other massive semi-detached or contact sys tems have been 



spect 


roscopically s 


tudied in both the UV and optical: 


LY Aur (HD 35921. Stickland et al. 


1994: 


Maver 


1968: 


Andersen. Batten. & Hilditch 


1974: 


PoDDer 


1982 


), 5 Ori A ( 


HD36486, 


Harvin et al. 


2002 


), V Puppis (HD 65818, 


Stickland et al. 1998 


Andersen et al. 


1983), and 


LZ Ccp (HD 209481. Howarth et al. 


1991: 


Harries. Hilditch. & Hill 


1998). All the Roche 



filling components have lower derived K values from the ultraviolet analysis when compared 
to the optical. We suggest that the difference between optical and UV semi-amplitude 



values for TU Mus originates in the close proximity of the two stars. In order to determine 
quantitatively how the tidal distortions affect the measured radial velocities, we need to 
complete the light curve analysis. From this analysis we will have a complete picture of the 
amount of distortion present. We save further discussion on the disagreement of optical and 
UV semi-amplitude values for §7 below. 



4. Projected Rotational Velocities 



We use a method developed previously (jPennyl Il996l ) to estimate the individual 
projected rotational velocities of the component stars from their cross- correlation functions 
(ccfs) with a narrow lined star. For HD 100213, we use HD 34078 as our template star. In 
the study above, we calibrated t he relationship between the ccf Gaussian width, and V sin i 



using the IConti &: Ebbetd (119771 ) data sample. The resulting calibration curve is given by 



V sin i = -3.830 x 10 -< V + 2.7903(x - 80.7, (1) 

The primary's and secondary's Gaussian widths of o = 164.4 ± 5.0 and 131.3 ± 5.0 km s _1 
correspond to a projected rotational velocities of V sini = 274.5 ±8.0 km s _1 and 219.6 ±8.0 
km s -1 . These agree within errors to those values determined by ST95 (250 ± 25 and 
195 ± 20 km s" 1 ) and AG75 (285 ± 29 and 240 ± 24 km s" 1 ). T03 do not quote errors on 
their Vsini values of 291 and 242 km s _1 , respectively. 
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5. Tomographic Reconstruction and Spectral Classification 



The separate primary and secondary star spectra can be derived using a tomography 
algorithm to reconstruct the individual spectra from the ense mble of composite sp ectra. 



Bagnuolo et al. 



1994). We 



19951 ). but 



The algorithm is described in detail in a separate paper (see 
have commented on some limitations of the algorithm previously (jThaller et al. 
we note here that strong wind features, such as the P Cygni lines, may be reconstructed 
incorrectly since their radial velocity curves can be very different from those associated with 
the stars themselves. Since the tomographic reconstruction is based on orbital Doppler 
shifts, the reconstruction will be ambiguous in the vicinity of such wind features. 

The tomographic reconstruction of the primary and secondary spectra is based on the 
assumed flux ratio which may not be well determined at the outset. The composite spectra 
are first separated using the ccf intensity ratios of r cc f = 0.49 as an initial estimate of the 
UV flux ratios. The resultant individual stellar spectra are then classified. For the spectra 
of the primary and secondary we use the spectral diagnostics listed below. We can then 
determine the true UV intensity ratio by numerical tests involving the cross-correlation of 
the template with a simulated binary spectrum formed from spectra of stars of the same 
spe ctral types as the p rimary and the secondary, with a range of input flux ratios (for details 



sec 



Penny et al. 



19971 ). For the case of TU Mus (07 V + 08 V), our measured ccf intensity 



ratio is 0.49 ± 0.05 using the template star, HD 34078 (09.5 V). We found that we could 
match this ratio by combining the spectra of the single stars HD 36879 (as the primary, 07 
V) and HD 48279 (as the secondary, 08 V) using a UV flux ratio, r uv = 0.48 ± 0.05. 

Our method for estimating the spectral types and luminosity classes of the separated 
primary and secondary spectra is b ased on the equivalent width m easurements of several 



UV absorption lines (for details see 



Penny. Gies. fc Bagnuolo 



19961 ). We identify the two 



stars as 07 V (06.5 - 07.5) and 08 V (07.5 - 09) based on an evaluation of the following 
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criteria: for the primary the equivalent widths of four lines (Si III A1299, He II A1640, Fe IV 
A1681, and Fe IV A1765); and for the secondary the same set plus Fe V A1429 and the line 



ratio of He II A1640/Fe V A1429. These classifications ar e, for the most 



with Walborn's (1971) composite classification of 08.5 V. 



j art, c onsistent 



Hilditch fe Bell 



(119871 ) estimated 



the types as 07.8 + 08.2, which a gain agree well with our determin ations. These spectral 



types correspond to temperatures (IMartins. Schaerer. fe Hillier 
kK and T cff s = 34.7 ± 1.5 kK. 



20021 ) of T effiP = 37.2 ±1.5 



6. Light Curve Analysis and Masses 



19971 ) shows two equally spaced minima of what 



The HIPPARCOS light curve fjESAI 

appears to be 0. m 54 0. m 48 in HIPPARCOS magnitude. We obtained, in electronic 

form, the AG75 differential photometric data from J. V. Clausen and J. Andersen. These 
observations show that in fact the primary eclipse is slightly deeper, 0.58 magnitude (in V). 
Both sets of data also have ellipsoidal variations present between eclipses due to the tidal 
distortion of the stars inside the Roche surfaces. We note that the time of minimum light 
for the HIPPARCOS observations is shifte d in phase by 0.02 8 compared to tha t of AG75. 



This is part of a larger trend examined by 



Maverl (12004h and 



Qian et al. 



(120071 ) suggesting 



that the period of TU Mus is lengthening due to slow mass transfer from the secondary to 
the primary. For the purpose of simultaneously fitting the HIPPARCOS and AG75 data, we 
calculated the phases of the photometric data using our final spectroscopic period and the 
time of primary minimum from the individual data sets themselves (i.e., the HIPPARCOS 
epoch of minimum light was taken as the time of the faintest observation). The differential 
V band photometry of AG75 was transformed to HIPPARCOS magnitude system (almost 
identical to V for hot stars) and are plotted in Figure 2. Both the shapes and depths of 
the eclipses and the ellipsoidal variations are dependent upon the stellar radii and orbital 
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inclination. 



19721 ) to 



We used the light curve synthesis code GENSYN ( IMochnacki fe Doughtyl 
produce V differential light curves. As T03 mention there is a faint visual companion to the 
system 15" to the SW which was included in both the HIP 'PARC 'OS and AG75 photometry. 
The visual magnitude of the companion from the Guide Star Catalog is V = 13.46. This 
represents 0.8% of the measured light, which we removed before fitting the light curve. 
The orbital parameters were taken from ou r above s olutio n, an d the physical fluxes and 



limb darkening coefficients were taken from iKuruczl (119941 ) and 



Wade fc Rucinskil (Il985l ) 



respectiy e 



y. W e neglect any treatment of radiation pressure following the suggestion of 



Howarthl (119971 ) . Rotation rates were determined using the model inclination and radii 
plus the projected rotational velocities given above. Making a fit of the observed light curve 
requires estimates of many parameters, and initially we decided to constrain all the orbital 
parameters according to our spectroscopic orbit, set the stellar temperatures according to 
the spectral classifications (T e fj p = 37.0 kK, T e ff s = 35.0 kK), and determine the ratio of 
radii from the visual flux ratio, F s /F p (determined from the method described in Penny et 
al. 1997 to be 0.49 ± 0.05). Then, each model light curve is a function of the remaining 
two parameters, inclination and primary polar radius. There are three characteristics of the 
light curve that we sought to match: eclipse depth, eclipse duration, and amplitude of the 
ellipsoidal variation. In order to fit the observed difference in depth between the primary 
and secondary eclipse, the primary's (secondary's) temperature was increased (decreased) 
to its highest (lowest) possible value, T e fj iP = 38.7 kK, T c s tS = 33.2 kK. The critical issue 
in matching eclipse depths is the temperature difference between the two stars, not their 
actual values. (Of course the actual values of T eS do have a large effect upon the calculated 
absolute magnitude of the system.) The ellipsoidal variations are a sensitive function of 
assumed stellar radii, we found that models with the correct ellipsoidal variation depth 
and eclipse duration and depth had inclinations 76?7 < i < 78?7. The best-fit model light 
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curve (for i = 77°. 7) is illustrated in Figure 2. In this model the stars are just in contact, 
both filling their Roche surfaces. The larger inclination value of 78?7 fit the eclipse depths 
very well, but the ellipsoidal variations were too small. Below 76?7 the widths of the 
eclipses become too broad and shallow. At the higher inclination, both stars are just inside 
their Roche volumes and at the lower, they are in slight overcontact. The presence of the 
ellipsoidal variations in conjunction with eclipses make this inclination determination very 
precise. 

EDITOR: PLACE FIGURE M HERE. 

Our best-fit inclination value agrees extremely well with the inclination presented by 
T03, who find an i = 77?8 ± 0?1. We caution against their small quoted errors. No doubt 
these are statistical in nature, however they can be misleading. To demonstrate we plot 
the secondary eclipse with our nominal best fit model (i = 77? 7), and our maximum and 
minimum inclination model fits (i ma x — 78?7, % m i n = 76?7) in Figure 3. All three fit the 
data very well. Owing to the smaller K values from our spectroscopic orbit, our estimates 
for the resulting separation and radii are quite different from that of T03. The radii at the 
best fit inclination are R p /R & = 7.2 ± 0.5 and R s /R & = 5.7 ± 0.5. Based on our V sini 
values, the equatorial rotation speeds of the two stars are 281 ± 15 km s" 1 and 225 ± 15 
km s _1 . The synchronous rotation rates, calculated using the equatorial radii, are 261 ± 8 
and 208 ± 7 km s -1 , for the primary and secondary, respectively. These agree within our 
errors, and both stars are rotating at or slightly above synchronous. To determine the 
distance to TU Mus we use the combined absolute magnitude of the pair from our light 
curve model, My = —4.6 ±0.1. We note that T03's use of absolute magnitudes from a 
calibration of spectral class to absolute magnitude is far less accurate. In fact their own 
light curve analysis produces log L/L Q values (4.8 ±0.2 and 4.5 ±0.2) that are incompatible 
with their adopted absolute magnitudes (which correspond to log L/L Q values of 5.3 and 
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4.9). We adopt a visual absorption value of l. m 2 and an apparent visual magntude of 8. m 23 
in accordance with AG75. We determine a distance of 2.1 ± 0.1 kpc, which agrees within 
errors with that from AG75, but which is much smaller than 4.8 kpc foun d by T03. O ur 



distance places the system in the nearby portion of the Carina spiral arm (JRusseil 



2003 ). 



EDITOR: PLACE FIGURE H HERE. 



7. Tidal and Irradiation Effects on Spectral Lines 

From our combined spectroscopic and photometric analysis we have a complete model 
of the TU Mus binary. We can now attempt to quantify why the optical and UV radial 
velocities result in disparate K values. We considered several effects. The first is the shift 
of the the stars' centers of light from their centers of mass. We used GENSYN to produce 
synthetic line profiles of typical UV and optical lines at photometric phase = 0.25, with 
a center of mass velocity of zero. At this orbital phase, the absolute values of the radial 
velocities of the component stars should be equal to their semi-amplitude values. To create 
model line profiles, GENSYN requires unbroadened line profiles for the both the primary 



and secondary. We adopt He I A4471 as our typical 



optical line. The 



synt 



retic input flux 



Lanz &: Hubenyl (120031 ) using effective 



profiles were obtained from the OSTAR2002 grid of 
temperatures from our best fit light curve model, T e g tP = 38700 K, T e fj jS = 33200 K, and 
logg = 4.0 for both stars. Our measured UV radial velocities are not from individual line 
profiles, but from a cross correlation of the spectrum of HD 34078 with the observed binary 
spectra. We cross-correlated UV spectra from OSTAR2002 with our template HD 34078 
and used the resulting ccfs as our input UV line profiles. Radial velocities were obtained by 
fitting double Gaussians to the inverted UV and optical GENSYN created profiles and are 
presented in Table 3. Neither the UV nor the optical appear to be significantly affected by 
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the non-sphericity of their shape (Case 1). The resultant velocities although slightly smaller 
are well within the measuring error of the expected semi-amplitude values. 

EDITOR: PLACE TABLE E HERE. 

Next we considered how the inner hemispheres of the components are heated by 
irradiation from the companion star, by setting the number of reflection iterations in 
GENSYN to equal five (Case 2). This had again a small effect on the measured radial 
velocities and not in a manner that addresses the K disparity, as both the UV and the 
optical lines are shifted towards the binary center of mass. This is not surprising given 
that at both wavelengths higher temperatures mean increased luminosity and a greater 
shift of the line center. Next we considered whether absorption lines formed primarily 
in hotter regions, such as UV lines, would more preferentially form in these inner hotter 
hemispheres, resulting in lower K values. An additional input to GENSYN does allow us 
to specify how the primary and secondary line profiles vary with increasing temperature. 
We again used the OSTAR2002 grid to quantify how the strength of the input lines change 
with temperature and we created an additional set of binary synthetic profiles with all 
three effects (Roche geometry, irradiation, line variability). A temperature increase of 1000 
degrees is expected to cause a decrease in He I A4471 equivalent width of 9% while the 
predicted change is a decrease of 3% for the UV ccfs. Disapointingly, the changes are too 
small to alter significantly the measured radial velocities. 

One last effect we consider is how the heating of the inner hemisphere reduces the 
temperature gradient in the upper atmosphere. Line depth is a function of the difference 
in temperature between the depths where the continuum and the line are formed, and in 
general, the larger the difference, the deeper the line. UV lines are high excitation features 
that are formed at depths comparable to the continuum. Optical lines, on the other 
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hand, are lower excitation features that form at much higher levels. But in cases where 
irradiation from above occurs, the temperature where optical lines form will be higher than 
normal, thus reducing the line contrast. In this situation, optical lines formed on the inner, 
illuminated hemisphere will appear weaker while lines formed on the outer hemisphere have 
normal strength, and this shifts the effective line center away from the center of mass, 
leading to larger K values. This effect would not so greatly change the UV line contrast 
as they are already formed at temperatures close to that of the continuum depth. We 
suggest that this may be the cause of the larger semi-amplitude values from optical lines. 
Without an accuate model of how irradiation changes the temperature gradient in the 
upper atmosphere it is difficul t to quantify this effect. Howev er in a previous study on 



the Struve-Sahade (SS) effect, 



Gies. Bagnuolo. fc Penny! (119971 ) investigated the effects of 



increased wind-blanketing on the temperature gradient, and subsequent line strengths, for 
the secondary star in AO Cas. Lines from transistions that tend to decrease in strength with 
temperature were not as strong, while those transitions that increased with temperature 
were either unchanged or intensified. While the UV and He II lines are somewhat affected 
by the increase in external heating, the largest impact is in the He I features. The largest 
effect can be seen in their Figure 5, where the He I A6678 feature is drastically diminished. 
In their case, they argue that a colliding bow shock region near the surface of secondary, 
creates heating on the approaching hemisphere of that star which is seen as the SS-effect. 
For our situation, only the inner regions of the stars would be heated from irradiation. 

We find observational support for this hypothesis in the recent paper by L07 in their 
measured equivalent width values for two He I lines, A4026 and A4471. In both cases the 
lines variations are phase locked, with maximum values at phases when we see the outer 
hemisphere of the stars. This ia consistent with the idea that the optical lines formed 
where irradiation from the companion is significant will be weakened by the decrease 
in atmospheric temperature gradient. L07 also suggest that there may be a similar but 
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opposite effect in the UV spectra. In order to determine if this is indeed the case we 
measured the 'equivalent widths' of our CCF primary and secondary peaks at phases out of 
eclipse, see Figure 4. If contributions to the UV ccfs were larger from the inner irradiated 
regions we would expect the primary to appear strongest at phase 0.25 and the secondary 
at phase 0.75. Unlike the He I lines there does not appear to be a correlation between line 
strength and orbital phase, which suggests that the UV lines form more or less uniformly 
around th stars. 

It is difficult to determine the actual optical depth of formation for any particular line, 
but they should scale with excitation and ionization potentials. The He I lines examined by 
L07 are from electron transistions from the first excited level (with an excitation potential 
of 21 eV) and these have the largest K values of 252 and 375 km s" 1 . The ionization 
potential between He I and He II is slightly higher at 24.6 eV and from the L07 analysis, 
the K values for He II lines are lower at 246 and 353 km s -1 . Finally, the UV spectra from 
which we obtain our semiamplitudes, 214.5 and 345.8 km s -1 , are predominantly populated 
with lines from Fe IV and Fe V with ionization potentials of 30.7 and 54.8 eV. It appears 
that the use of low excitation lines in contact systems where irradiation effects are present 
leads to overestimates of the semi-amplitude values and hence to the overall size and masses 
of the system. It is clear from this analysis that the UV derived orbital elements are more 
representative of the true dynamic motion of the system. 

EDITOR: PLACE FIGURE ?? HERE. 

8. Placement in the H-R diagram and Discussion 

Our goal in this research has been to compare the masses obtained from observations 
with those found from theoretical evolutionary tracks, however since both stars of TU Mus 
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are Roche filling, we should not necessarily expect agreement with single star tracks. In 
this case we will also compare the component stars to models produced from interacting 
binary models. The observational masses can be determined from the combination of results 
from the spectroscopic and photometric orbits (m sin 3 ? and i). Masses from evolutionary 
theories require the individual temperatures and luminosities of the two stars, i.e., placing 
the individual component stars in the H-R diagram. We can estimate temperatures from 
calibrations of spectral type, and we can calculate luminosities from the observed flux ratio 
and the absolute magnitude of the binary (from the stellar radii associated with fits of the 
light curve). Individual visual magnitudes of the component stars can be obtained from 
the visual flux ratio of the two stars and the adopted My above. Once individual visual 
magnitudes are obtained, we con vert them to bolometric magnitudes and luminosities 



with bolometric corrections from 



Howarth fc Prinjal (Il989l ). Our calculated visual flux 



ratio is ry = 0.49 ± 0.05. The luminosities and radii derived are log L p /L Q = 5.02 ± 0.08, 
log L S /L Q = 4.55 ±0.09, R P /R Q = 7.2 ±0.5, and R S /R Q = 5.7 ±0.5. Our results are plotted 
in th e H-R diagram in Fig ure 5 together with evolutionary tracks for single massive stars 



from 



Schaller et al. 



(Il992l ). At this point in our analysis we typically compare the stellar 



positions with the evolutionary tracks and determine what masses they would predict. 



EDITOR: PLACE FIGURE ?? HERE. 



Combining the light curve model with the results from the spectroscopic orbit results 
in observational masses of M P /M Q = 16.8 ± 0.4 and M s /M & = 10.5 ± 0.3. Masses predicted 
from single star evolutionary tracks are significantly higher, 27.7±1.5M and 18.O±1.5M . 
Of course since both stars are in Roche lobe contact, we would not expect them to have 
evolved as single objects. The severe overluminosity of these Roche filling stars is in 
agreement with other interacting massive binaries, specifically HD 115071 , HD 36486 , 
HD 35652, HD 209481, BD+66°1521, HD 228854, HD 190967, HD 106871, and HD 35921. 



These systems are listed in Table 4, along with their individual periods, mass ratios, 
component spectral classes, evolutionary masses from single star tracks, and dynamical 
masses. 



EDITOR: PLACE TABLE g| HERE. 



An interacting massive binary is a complicated beast. Trying to determine how it got 
to its current state is a little like trying to unspill milk. One method is to compare its 
current state with other similar systems. Interacting massive binaries fall into two broad 
categories: semi-detached and contact. Its interesting to note that in all but one of the 
semi-detached systems in Table 4, it is the secondary that is in RLOF. In the Case A (donor 
star still in core hydrogen burning, CHB) scenario, the radius of the initially more massive 
star slowly reaches its Roche surface as it evolves along the main sequence. The initially 
lower mass star, due to its slower evolution, is not the first to overflow. Applying this to 
the semi-detached systems we infer that the current mass loser (or secondary) was initially 
the more massive object. The period of a binary will shrink during mass transfer as long as 
the more massive star is losing mass. After enough mass has been transferred to 'flip' the 
mass ratio, the separation between the stars will increase allowing RLOF to cease. That we 
mainly se e systems with the less massive object in RLO F agrees with theoretical predictions 



(WLB01, 



Wellstein. Langer. fc Braun 



2001 



Pols 



19941 ) that the initial transfer of mass 



is exceedingly quick (fast Case A), but that the turning off of slow mass transfer (slow 
Case A) may take a considerable amount of time. The mass ratios of these semi-detached 
systems range from 0.353 — 0.556, indicating that there has been significant mass lost 
from the intially more massi ve star. Unfort unately only one of the semi-detached systems, 
HD 209481, was observed by [Thaller (199?]) in her survey of Ha emission in close, O-type 



star binaries. It did not have strong the Ha emission that is associated with large RLOF 
gas streams. Again this suggests that although the secondary is filling its Roche volume, 
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there is not a significant amount of mass being transferred at this time. 

If the semi-detached systems are the expected result when the intially more massive 
star in a binary evolves to reach its Roche radius, what scenario leads to a contact system? 
The are some obvious differences between the two groups. The observed mass ratios for 
the three contact systems in Table 4 are slightly higher (q av = 0.64) than those of the 
semi-detached systems (q av = 0.47). The components are also more closely matched, within 
one spectral type of each other and with similar luminosity classifications. WLB01 present 
evolutionary calculations for 74 close binary systems, with masses ranging from 26...6M 
for Case A and B scenarios. On the basis of our derived luminosity classes (both stars are 
V) and the log g values determined from our light curve analysis (3.98 and 3.97), we expect 
that HD 100213 is undergoing Case A evolution. According to the WLB01 models, Case 
A contact systems can occur when the initial period of the system is very small and/or 
when the mass ratio, q = M 2 /Mi, is either very small or very large. For a short period 
system, the initially more massive star reaches its Roche radius fairly early in its evolution. 
If its companion has much less mass than itself (q << 1) the ensuing mass transfer will 
proceed very rapidly as the exhange of any mass reduces the period causing the loser star 
to overflow at a faster rate. The combination of the large increase in mass of the gainer star 
and the smaller orbit causes it to expand and reach its Roche volume. If the companion has 
a mass similar to the loser star (g « 1), it will gain a much smaller amount of material from 
mass transfer. In this case only a small mass transfer will 'flip' the mass ratio. However the 
increase in mass will cause the gainer star to 'rejuven ate' . Rejuvenation is when the star 



adap ts its convective core mass to its new total mass (IHellings 



1983 . 



1984 



Braun fc Langer 



19951 ). The new larger mass gainer star begins to evolve at a faster pace than its companion 
and reaches its Roche limit. The state of the gainer star (original secondary) when it reaches 
its Roche volume determines whether the contact state is called Case AA (secondary still 
in CHB) or Case AB (secondary between CHB and CHeB). In Case AB the mass transfer 
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rate is very large as the star is rapidly increasing in radii. This exchange happens quickly 
and will increase the orbital period and separation significantly 

For HD 100213 we suggest the following scenario. The two stars began with a mass 
ratio very close to unity and a period slightly smaller than the current value. The actual 
initial masses and period depend upon how much of the lost mass is transferred. The 
initially more massive star reaches its Roche volume and mass is transferred to the gainer 
star. The stars move closer together for a brief time until the gainer star is more massive. 
Mass transfered after this point causes the stars to move further from one another, however 
this does not cease the RLOF. Fast Case A mass transfer continues and doesn't end until 
the loser star is much less massive. WLB01 explain that since the M core /M enve i ope is 
increasing for the loser star this causes its radius to increase. This agrees well with our 
observed mass ratio of 0.625. Finally the gainer star expands due to the increase in mass 
and its rejuvenated interior until it also reaches its Roche v olume. This is current status 



of the binary. W hat is the curren t state of mass tra nsfer? 



historical records (jOosterhofi 



1928 



1930 



Knipe 



1971 



Mayer 



(120041 ) examined the 



AG75, ST95, T03) and argue t hat the 



Qian et al. 



)7) 



(20071) found 



period of TU Mus lengthens, but not regularly. Another study by 
a long term increase in period along with superimposed cyclical variations. The cyclical 
variations they attribute to the presence of a nearby (~ 50 AU), low mass (~ 2M Q ) tertiary 
companion. They argue that the long term increase in the observed periods is consistent 
with a conservat ive mass trans fer rate of 4.2 x 1O" 7 M yr _1 from the lower mass secondary 



to the primary. iThallerl (jl997l ) observed HD 100213 looking for Ha emission, which we 
associate with mass transfer, and detected none. We suggest that the system is just barely 
in contact, with both stars still in a CHB stage, Case AA. 

We highly recommend that this rare contact system, now with well determined mass es, 



radii, separation, and temperatures be modelled using binary evolutionary codes. 



Pols 
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( 119941 ) also examined Case A and B interacting systems and determined that those with 



small initial periods and small or lar ge m a ss rat ios will reach contact. However due to 



his differing treatment of convection, 



Polsl (119941 ) predicts that all gainer stars will reach 



contact through Case AB. All gainer stars are quickly rejuvenated to the point that they 
reach the end of CHB and then expand rapidly to contact. Another issue this system can 



address is how much mass lost from the loser is accreted upon the gainer. Both 



Polsl (I1994T ) 



and WLB01 have treated the mass transfer as conservative, i.e., all mass lost by the donor 



is accreted by the gainer. The mechanism for non-conservative mass trans: 



spinning up of the gainer star unti 



Langer. Wellstein. fc Petrovic 



2003 



it reaches its critical rotation velocity ( Wellstein 



Langer et al. 



2004 



Petrovic fc Langer 



'er is the rapid 



2001 



20041 ). Neither 



star in this system i s nea r critic al rotation, which suggests that the mass transfer has been 



conservative. What iPolsl (Jl994j) and WLB01 agree upon is the prediction that Case A 
contact systems will end in a binary merger. Both find that the resulting periods are too 
small for successful common envelope ejection. A merger results in a very exotic object, 



which may be a precursor to Gamma Ray Burster or other strange supernova. 
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Figure Captions 

Fig. 1. — Radial velocity measurements for the primary (filled circles) and secondary (open 
circles) based on IUE spectra. Solid lines are drawn for both the primary and secondary 
radial velocity solutions. 

Fig. 2.— The observed HIPPARCOS light curve (filled circles) and AG75 differential V 
light curve (filled triangles offset by +5.74 mag) plotted together with the optimum model 
fit (after adjusting the stellar effective temperatures, solid line). 

Fig. 3.— The observed HIPPARCOS and AG75 differential V light curves (offset by +5.74) 
plotted with our minimum acceptable inclination model (short dashes), maximum acceptable 
inclination model (long dashes) and optimal model (solid line). Only the secondary eclipse 
is shown. The data points are greatly reduced in size to better see the model fits. 

Fig. 4. — Measured 'equivalent widths' of our CCFs outside of eclipses for the primary (filled 
circles) and secondary (open circles). 

Fig. 5. — A HR diagram of the binary system TU Mus. The filled (open) circle represents 
the primary (secondary) component. The overdrawn evolutionary tracks are for single stars 
from Schaller et al. (1992). 



Table 1. IUE Radial Velocity Measurements 



SWP 


HJD 


Orbital 




(0- 


~C)i 




v 2 


(0 - -C) 2 




Image 


(-2,400,000) 


Phase 


(km s _1 ) 


(km 




W 1 


(km s _1 ) 


(km s _1 ) 


w 2 


15916 


44970.915 


0.134 


154.4 




20.5 





-214.3 


11.6 





19642 


45429.986 


0.048 


199.5 




3.9 


1 


-330.0 


-4.6 


1 


19643 


45430.049 


0.093 


183.4 




13.9 


1 


-274.3 


9.0 


1 


19651 


45430.480 


0.404 


-207.5 




-21.7 


1 


284.8 


-4.7 


1 


37879 


47883.732 


0.792 


92.0 




45.6 





-119.2 


-34.3 





37880 


47883.809 


0.847 


147.8 




33.9 





-241.1 


-47.5 





37906 


47887.590 


0.573 


-194.3 




7.5 


1 


313.0 


-2.2 


1 


37907 


47887.654 


0.619 


-171.9 




-5.3 


1 


282.4 


23.8 


1 


54353 


49817.760 


0.906 


165.3 




-3.2 


1 


-302.3 


-20.7 


1 


54356 


49817.933 


0.030 


191.3 




-10.2 


1 


-325.7 


9.0 


1 


54360 


49818.117 


0.163 


130.6 




28.1 





-234.6 


-59.3 





54371 


49818.514 


0.449 


-206.8 




6.0 


1 


311.9 


-21.2 


1 


54374 


49818.700 


0.583 


-189.1 




6.2 


1 


316.2 


11.4 


1 


54378 


49818.917 


0.739 


-64.1 




-40.5 





64.6 


36.4 





54382 
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0.869 
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16.0 
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-54.2 





54390 


49819.481 
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22.6 





-244.4 


-39.0 





54394 


49819.667 
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-126.3 




-76.7 





83.1 


13.3 





54397 


49819.813 


0.385 


-194.7 




-24.3 





274.7 


10.0 





54407 


49820.633 


0.976 


200.0 




-2.9 


1 


-332.3 


4.7 


1 


54411 


49820.805 


0.100 


167.5 




3.5 





-276.5 


-2.1 





54415 


49820.994 


0.237 


42.1 




33.3 





-87.4 


-63.2 





54419 


49821.185 


0.374 


-181.9 




-21.6 





271.0 


22.7 
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Table 1 — Continued 

SWP HJD Orbital V 1 (O - -C) 1 V 2 (O - -C) 2 

Image (-2,400,000) Phase (km s^ 1 ) (km s" 1 ) W 1 (km s" 1 ) (km s" 1 ) W 2 

54426 49821.471 0.580 -191.7 5.3 1 303.9 -3.6 1 



-29- 



Table 2. Circular Orbital Elements 



Element 



ST95 



T03 



Primary only Secondary only 



P (days) 



1.3872827(17) 1.38728653(2) 



T (HJD 2,400,000+) 49817.879(3) 49817.8820(1) 



K 1 (km s- 1 ) 
K 2 (km s- 1 ) 

V (km s- 1 ) 
V i (km s- 1 ) 
V 2 (km s- 1 ) 

mi sin 3 i (M ) 
m 2 sin 3 ? (M ) 

ai sin i (i? ) 
a 2 sin i (i? ) 

r.m.s. (km s _1 ) 
r.m.s.i (km s _1 ) 
r.m.s.2 (km s _1 ) 



216.7(2.7) 
345.4(3.1) 

-12.7(2.1) 
-11.8(2.1) 
-14.2(3.6) 

15.70(34) 
9.85(24) 

5.94(7) 
9.47(8) 

13.0 



249.0(2.8) 
382.6(4.3) 

-4(4) 



21.9(8) 
14.3(4) 

6.97(8) 
10.5(1) 



1.387282(7) 
49817.662(5) 
214.5(2.4) 



-9.3(2.1) 



15.64(44) 



5.88(7) 



10.1 



1.387280(6) 
49817.666(5) 

345.8(3.2) 



4.9(2.9) 



9.71(32) 



9.48(9) 



13.7 



-30- 



Table 3. V r from Synthetic Profiles 



UV K, P UV V ryS Optical V rjt Optical V Ty 



Case 


Reflection Effects? 


Line Variability? 


(km s l ) 


(km s l ) 


(km s v ) 


(km s 


Geometric 






-214.5 


345.8 


-214.5 


345.8 


f 


off 


off 


-210.9 


344.9 


-209.9 


345.1 


2 


on 


off 


-210.6 


343.1 


-209.5 


343.7 


3 


on 


on 


-211.6 


343.3 


-209.7 


345.2 



Table 4. Interacting OB Star Binaries 



Star 


P 


1 


Primary 


Secondary 




(M ) 




(M ) 


Roche Filling 




Name 


(days) 


(Af 2 /Mi) 


Type 


Type 


Primary 


Secondary 


Primary 


Secondary 


Star? 


Reference 


Scmi-Dctachcd Systems 


HD 35652 = IU Aur 


1.81 


0.503 


B0.5 V 


B0.5 V 


15 


13 


14.5 


7.3 


Secondary 


2 


HD 115071 = V961 Cen 


2.73 


0.575 


09.5 V 


B0.2 III 


18 


15 


11.6 


6.7 


Secondary 


4 


HD 209481 = LZ Ccp 


3.07 


0.415 


08.5 III? 


O9.5-B0 V? 


21 


18 


15 


6.3 


Secondary 


2 


HD 106871 = AB Cru 


3.41 


0.353 


08 V 


B0.5 ? 


31 


16 


19.8 


7.0 


Secondary 


3 


BD +66° 1521 = XZ Cep 


5.10 


0.405 


09.5 V 


Bl III 


15 


14 


15.8 


6.4 


Secondary 


1 


HD 36486 = <5 Ori A 


5.73 


0.500 


09.5 II 


B0.5 III 


28 


12 


11.2 


5.6 


Primary? 


5 


HD 190967 = V448 Cyg 


6.52 


0.556 


09.5 V 


Bl Il-Ia 


15 


14 


25.2 


14 


Secondary 


1 










Contact Systems 












HD 100213 = TU Mus 


1.38 


0.652 


07 V 


08 V 


28 


18 


16.8 


10.5 


Both 


10 


HD 228854 = V382 Cyg 


1.89 


0.742 


07.3 V 


07.7 V 


32 


24 


26 


19.3 


Both 


1 


HD 35921 = LY Aur 


4.00 


0.528 


09 III 


09.5 III 


34 


28 


24 


12.7 


Both 


6,7,8,9 
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